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C| Abstract 

In the civilian aviation field, the radar detection of hazardous weather 
phenomena (winds) is very important. This detection will allow the 
avoidance of these phenomena and consequently will enhance the safety 
of flights. In this work, we have used the wavelets method to estimate 
the mean velocity of winds. The results showed that the application 
of this method is promising compared with the classical estimators 
(pulse-pair, Fourier). 
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1 Introduction 

Mainly in civil aviation, the meteorology plays a big role in the security of 
nights. Aeronautics depends on reliable information related to the present 
and future meteorological conditions. This information is given by a me- 
teorological radar. This last detects and locates even distant atmospheric 
disturbances in form of signals in microwaves. 
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Several methods of signal processing are worked out, others offered, to 
complete successfully the extraction of useful information to be transmitted 
to the user. 

The pulse-pair method acts on the Doppler radar signal in the temporal 
domain by autocorrelation. As for the wavelets method, it spreads out in 
the spectral domain while having references to time. 

The main purpose of this work is to improve the estimates of the average 
velocity and variances of the detected meteorological disturbances. The use 
of an algorithm of denoising based on the wavelet transform will lead to this 
improvement. 

In section 2, we will introduce the estimation theory of spectral moments. 
This estimation will be made by two estimators, one temporal, named pulse- 
pair and another one scalo-temporal, the wavelets. The section 3 introduces 
results and comments linked to calculations and simulation. Finally a gen- 
eral conclusion will be presented in section 4. 

2 Estimation theory 

The first three (03) moments of the power spectral density of the Doppler 
spectrum are directly linked up with the desired atmospheric basic param- 
eters. The radar reflectivity (Z), the radial velocity (Vr) and the spectral 
width of velocities (W). [1], [2], [3], [8]. 

Since the return signal in the radar from a range cell (space dimension of 
an impulse) is generated by the back scattering of a big number of randomly 
distributed particles and/or by variations of the refraction index of the at- 
mosphere/air, then the process of the received signal is considered (central 
limit theorem) or approximated by a gaussian random process. 

The shape of the received signal reflected by all particles met by the 
transmitted beam is: [3], [5], [8]. 

V(t,R) = Y,^ f ^- 2 ^) .w{t - 2^) (1) 

i 

with: 

Af. Complex amplitude of the individual particle i, 
ff Doppler frequency of the particle i, 
Ri'. Distance from the radar to the particle i, 
w[t — 2—): Weighting Function of distance. 

At the k th To the kth radar impulse, the received wave shape may be 
written in the complex form [3], [5], [6] 
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V k (R) = I k (R) + jQ k (R) 



(2) 



with: 

Ik(R)' In phase component. 
Qk(R)' Quadrature component. 

The estimation of the meteorological parameters is typically accom- 
plished in the signal processing on the basis of range cells. 

2.1 Pulse-Pair Method 

The working frequencies of a meteorological radar are in the order of 10 9 Hz 
corresponding to wavelength of some centimeters (jlO cm). Generally, the 
meteorological targets move with speeds lower than 50 m/secs. The Doppler 
Effect would translate these speeds into a shift of the transmitted frequency. 
This shift is a few hundred of Hz: much too weak to be directly measured. 
To reach it, we measure the phase shift between the return of two successive 
impulses having probed the same volume of space. Indeed, after the return 
of the second impulse, the target would have changed position which would 
be translated in a phase shift between both impulses. Hence the pulse-pair 
method is derived. 

2.2 The pulse-pair technique 

The classical pulse-pair technique gives the first two spectral moments (speed 
Doppler and its variance) of the weighted Doppler spectrum of the reflectiv- 
ity from the function of autocorrelation of the radar complex signal [1], [2], 



Components I and Q are supposed to be independent statistically [5, [6], 



The function of autocorrelation can be expressed for the sequence V[n] 
of length N and whose samples are spaced by Ts, by: [3], [5], [6] 



[3], [4], [5]. 



V[n] = I[n] + jQ[n] 



(3) 



[7]- 




N-2 




n=0 



The mean velocity of a weather target is estimated by [5] : 



A 



arg(R[l\) 



(5) 



v = — 



4ttT s 
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And consequently, the spectral width of speeds is done by: 

A 



w 



2y/2irTs 



'In 



Rs+n[0] -R n [0] 



R[l] 



(6) 



Where index S+N indicates the weather signal merging in the noise N. 



3 Wavelet Method 

The fundamental assumption in signal 
echo signal consisting of two parts. 



processing is to consider the radar 



y = fa + e (7) 

Where f a produced by the Gaussian process of atmospheric diffusion and 
e a noise resulting from various but mainly thermal sources [1]. 

In the following, we will deal with the problem by using the wavelets to 
filter the collected signal then apply the already seen method of processing 
(pulse-pair) to extract the various spectral moments of interest. 

The primary reasons of the effectiveness of processing by wavelets are 
linked to the nature of the noise: 

- it is nonstationary 

- it is transitory and of unknown form 

We will suppose that the Gaussian model fits the atmospheric signal f a 
and noise e. 



3.1 Discrete wavelet transform 

The dyadic discrete wavelet transform (DWT) of a signal y(t) is [9], [10] 

/+oo 
y(t)2~^^*(2~ m t-n)dt (8) 
-oo 

Where: 

m: represents the considered scale. 
n: represents the translation of the wavelet ^(t). 
T my n'- represents a coefficient of correlation between y(t) and *&(t). 
^m,n = 2~^ty*(2~ m t — n), represents the orthonormal basis of the 
mother wavelets. 

Consequently the reconstruction of y(t) from the coefficients T m ^ n is then: 
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3.2 Multiresolution analysis 

If one samples a continuous signal u with a regular interval unit, one will 
obtain a sequence with discrete values u n =u (n) , n € Z. [11] 

The dwt on one level breaks up the sequence u n =u (n) in two sequences 
al and dl by a low-pass filter h and a high-pass filter g, both followed by 
an undersampling of order 2. 

The sequence al obtained is called approximation and contains infor- 
mation of low frequencies of u n =u (n) , while the sequence dl contains in- 
formation of high frequencies of u n =u (n) called details, h and g are the 
decomposition filters of finite lengths (FIR) . 

a 1 = (u*h)±2 (10) 

d 1 = (u*g)l2 (11) 

The operation of undersampling by 2 means to take every other sample. 
It is very useful during the reconstruction operation. This one is carried out 
by means of filters of reconstruction h and g. The filters of analysis and 
those of synthesis are in quadrature mirror (QMF). 

Indeed, the sequence u n can be found by the expression: 

u n = (a 1 |2) *h+ (d 1 t 2) *g (12) 

The oversampling of order 2 inserts zeros in the sequences in order to 
find the initial sequence u n with the same initial number of samples. 

If the operation of decomposition is repeated J times, one obtains the 
algorithm of filter bank applied to the approximation outputs (see figure 1). 

In the opposite direction, the reconstruction of u n is carried out starting 
from the approximation a 4 and the details d 1 , d 2 , d 3 and d 4 . 

3.3 Thresholding 

The wavelet coefficients correspond to the details of a signal. A detail can 
be been ignored without affecting significantly the data. 

There exist various methods of thresholding of the wavelet coefficients: 
hard thresholding and soft thresholding. 

Hard thresholding: 
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Figure 1: DWT of « n to level J = 4: the sequence u n is decomposed into 
an approximation a 4 and details d 1 , d 2 , d 3 et d A 



d s = \* 31 H<M (13) 
s \ d \d\ > A J v ' 

Soft thresholding or shrinking specific to the complex signals: [12] 

*-<(i-£) (") 

As for the various thresholds, one can use: 

- The universal threshold: = v / 21og(A f /o5 (N)) 

- The minimax threshold: (see table I) 

- The SURE threshold: (Stein's unbiased risk estimate). 

Where d, A and N are the wavelet coefficients, the threshold and the 
length of the signal, respectively. 
Table I: Extract of table 1 of [12] 



N 


Minimax 

(An*) 


16 


1.763 


128 


1.973 


256 


2.176 


512 


2.371 


1024 


2.560 


2048 


2.741 



3.4 Signal processing algorithm 

The thresholding is efficient means to ignore the weakest details, compared 
to the selected threshold , A, which one can compare to the noise and pre- 



serve only the most important wavelet coefficients. One rebuilds, then, the 
signal from the remaining coefficients which represent the denoised data (see 
figure 2). 
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Figure 2: Synoptic diagram of the wavelet denoising method. 



4 Results and comments 

In this work, we have generated data representing the echoes of a weather 
radar of wind disturbance. 

For that, we have used the algorithm of ZrniA§ developed in 1975 and 
included in works of R.D. Palmer. These data consist of signals I and Q. 
These calculations were carried out on a set of ten range cells taken on one 
radar azimuth. 

In figure 3, a, b, we have represented signals I and Q and their spectra 
respectively for the range cells 1 and 5. As one can see it on the figure, this 
spectrum is a Gaussian form. 

The estimates of the mean velocity and the spectral width of the weather 
disturbances already generated above are represented on figures 4 and 6. 
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Figure 3: signals I and Q and Doppler spectrum of the cell Nl 

It is noticed that the estimate of the mean velocity of the wind by the 
pulse-pair algorithm is close to the real speed. On the other hand, the 
estimate by the method of Fourier is less close. 

We also notice that the estimation of the spectral width is less important 
with the wavelets than with the pulse-pair algorithm. 

In addition, by applying the wavelet method, we have obtained better 
results compared to Fourier and pulse-pair, see figure 4. This indicates the 
interest to have used this method. These insufficiencies of the algorithm of 
Fourier are inherent in his approach and they are due to: 

- Resetting the frequencies which are apart from the field of work, which 
generates spectral losses. 

- Spectral resolution (close frequencies). 

5 Conclusion 

In this work, we have re-examined the estimate of the spectral moments of 
the signals received by a Doppler radar, namely the moments of orders one 
(mean velocity of the wind) and two (spectral width). 

It comes out from it, through the results obtained, that the algorithm of 
denoising by the means of the wavelet transform of the signal to be treated is 
efficient. Our results were consolidated by the comparison of the errors of the 
various used algorithms: pulse-pair and wavelet. The proof is the reduction 
of the relative error made on the parameters estimated compared to the 
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Mean velocity with P.P and wavelets- SNR = 12 
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Figure 4: Estimate of mean velocities of the wind by the Pulse-pair and 
wavelet estimators. 



input data compared to the results obtained by the pulse-pair algorithm as 
well as the reduction of the spectral width 
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Spectral width by P.P and wavelets : SNR = 12 
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Figure 6: Spectral width obtained the Pulse-pair and wavelet estimators. 
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